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ABSTRACT 

High resolution spectra are necessary to distinguish and correctly measure the Balmer emis- 
sion lines due to the presence of strong metal and Balmer absorption features in the stellar 
continuum. This accurate measurement is necessary for use in emission line diagnostics, such 
as the Balmer decrement (i.e. HafHfi), used to determine the attenuation of galaxies. Yet at 
high redshifts obtaining such spectra becomes costly. Balmer emission line equivalent widths 
are much easier to measure, requiring only low resolution spectra or even simple narrow band 
filters and therefore shorter observation times. However a correction for the stellar contin- 
uum is still needed for this equivalent width Balmer decrement. We present here a statistical 
analysis of the Sloan Digital Sky Survey Data Release 7 emission line galaxy sample, using 
the spectrally determined Balmer emission line fluxes and equivalent widths. Using the large 
numbers of galaxies available in the SDSS catalogue, we determined an equivalent width 
Balmer decrement including a statistically-based correction for the stellar continuum. Based 
on this formula, the attenuation of galaxies can now be obtained from low spectral resolu- 
tion observations. In addition, this investigation also revealed an error in the H/3 line fluxes, 
within the SDSS DR7 MPA/JHU catalogue, with the equivalent widths underestimated by 
average ~0.35A in the emission line galaxy sample. This error means that Balmer decrement 
determined attenuations are overestimated by a systematic 0.1 magnitudes in Ay, and future 
analyses of this sample need to include this correction. 
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1 INTRODUCTION 

The Balmer lines are the most well known and observed emission 
lines in astronomy, being both strong lines in the optical and ubiq- 
uitous as they arise from recombination to the n = 2 level of the 
most common element, hydrogen. As the atomic structure of hydro- 
gen is so well understood, the strength of the emission line fluxes 
can be well determined if the radiation field ionizing the hydrogen 
gas is known, and, importantly, the relative fluxes of the resulting 
hydrogen emission lines are only weakly dependent on the local 
conditions. Given that the ratios are well determined, the difference 
between the measured ratios of the Balmer lines and the intrinsic 
values expected can be used to determine the reddening of galaxies, 
or, more accurately, the ionized regions within them. In association 
with a attenuation/reddening-law and selective-to-total attenuation, 
R v , the reddening can then give the total attenuation of a galaxy 
(e.g. the oft used work of Calzetti 2001). This possibility of us- 
ing the hydrogen emission lines, in particular the ratio of the two 
strongest Balmer lines Ha and H/3, to measure the reddening and at- 
tenuation has been known and utilized for many years (e.g. lBermanl 
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1 1936L was one of the first mentions of the Balmer decrement being 
affected by the path length through an absorbing medium). 

However, the accurate measurement of the Balmer decrement 
(i.e. the ratio of Ha/H/3) requires the measurement of both the rel- 
atively weak H/3 line an d the continuum underneath it to distin- 
guish the line. Separation of the underlying continuum from the 
Balmer emission lines is vital as the metal absorption lines and, es- 
pecially, Balmer absorption lines present in the spectrum of later- 
type stars act to weaken or hide the relative flux of the emission 
lines in the integ rated spectra of galaxies. Such was shown by 
iLiang et al] J2004t) . where the H/3 line was not observed in a se- 
lection of galaxies in low-resolution spectra (R = 150) from the 
the Canada-France-Redshift Survey, but moderate resolution spec- 
tra (R > 600) revealed the weak H/3 lines hidden by both dust 
and absorption lines. The exposure times needed to obtain suffi- 
cient spectral resolution to distinguish the emission lines from the 
underlying stellar continuum of galaxies means that studies of the 
Balmer decrement tend to be biased to high emission-line equiva- 
lent width objects. While this bias may not be a serious issue, with 
dustier, more attenuated objects tending to have higher specific star 
for mation rates and thus h igher emission line equivalent widths (see 
e.g. lda Cunha et alj20ld) . such biases do tend to limit samples for 
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investigations of dust and star formation. This is especially so at 
higher redshifts where the high spectral resolution needed to re- 
solve both the line and continuum limits surveys to the brightest 
objects. 

Even with moderate resolution spectra, the Balmer emission 
line fluxes are still sensitive to the way the stellar absorption is 
accounted for, and this can be a substantial source of error, particu- 
larly for weak emission line sources. 

It is these issues that motivates us to examine the possibil- 
ity of determining the decrement from the equivalent widths of the 
Balmer emission lines. Emission line equivalent widths do not re- 
quire high resolution spectra to distinguish the line from the contin- 
uum, allowing the use of low resolution spectra, such as at R~ 300, 
enabling fainter objects or more objects to be observed for the same 
exposure time as high resolution spectra. Even at only R~ 100, as 
will be available with the multi-object NIRSpec on JWST, the in- 
fluence of H/3 on the interpretation of galaxy spectra can still be 
significant (Pacifici et al., in prep), and the line is still detectable at 
a signal-to-noise ~ 5 for some of the brighter emission-line galax- 
ies. The issue at these low resolutions (i.e. R < 300) becomes one 
of distinguishing individual lines, such as [Nn]/i6584A from Ha, 
not detecting the lines. 

To examine how the Balmer line equivalent widths can be used 
to determine t he Balmer decrement we use the Sloan Digital Sky 
Survey (SDSS lAbazaiian et alj|2009l) . which contains a large spec- 
troscopic sample of emission line galaxies covering a wide range 
of galaxy types and properties, including attenuations. With such 
a wide range of galaxies, and reasonably high resolution spectra 
(R ~ 1900), the SDSS provides the perfect sample for examining 
this issue. With the large number of SDSS galaxies it is possible to 
obtain a statistically representative estimate for the correction fac- 
tors needed (i.e. the Balmer absorption lines) and thus determine 
the Balmer decrement from equivalent widths alone. 

This method should be seen as complementary to the stel- 
lar spectral synthesis continuum fitting methods widely used 
in the analysis of g alaxy spectra. Codes such as PLATEFIT 
dTremonti et al.l 1 20041 used for the determination of the stel- 
lar properties in the S DSS MPA/JHU catalogue used here) and 
(Lamareille et al. 2006, applied this to low er resolution, higher red- 
shift data in the VVD S sample), pPXF (jCappellari & Emsellerr] 
|2004|). STARLIGH T dCid Fernatides et al.l l2005h. S TECKMAP 
dOcvirk et al]|2006l) , and VESPA dToieiro et alj|2007l) . all use lin- 
ear c ombinations of synthetic stellar population spectra (such as 
from iBruzual & Charlotl 120031) to fit the full observed spectra of 
galaxies using various optimized maximum likelihood approaches. 
These codes have been created to extract the maximum possible 
information from gala xy spectra given de generacies and noise (see 
e.g. the discussion in lOcvirk et al.ll2006l) . and thus are clearly the 
best approach when strong continuum is detected. Yet the amount 
of possible information to be extracted reduces with both decreas- 
ing signal-to-noise ratio and spectral resolution. In addition, these 
methods are limited by the available spectral libraries, which may 
not cover the full parameter range needed to match the observed 
galaxies, and may have intrinsic issues, as we demonstrate here in 
an issue we discovered in the course of this paper. Thus an em- 
pirical method as we explore here is fully complementary to the 
spectral synthesis methods used in most works. 

In the following sections we introduce the Balmer lines, both 
in absorption and emission (^2), the SDSS emission-line galaxy 
sample(^3] provide a possible way to determine the Balmer decre- 
ment from the equivalent widths(©, and finally also point out an 
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Table 1. Balmer lines, including rest-frame wavelengths (Air), and their 
ratio s relative to H/3 for n e = 10 2 c m -3 and 3 different te mperatures ([values 
from Dopita & Sutherland] b003l) , based on data from IStorev & HummeJ 
1 1995)). 



interesting problem with the fitting of the stellar continuum in the 
SDSS (©. 



2 THE BALMER LINES AND DECREMENT 
2.1 Balmer emission lines 

The Balmer emission lines in the interstellar medium arise predom- 
inantly from the recombination and subsequent cascade of elec- 
trons to the n = 2 level of hydrogen. While collisional excitation 
can also contribute to t he Balmer line emission in hot media (see 
e.g. lFerlandetalj 2009). photoionization and recombination are the 
predominant energetic processes in most galaxies. 

As the atomic structure of hydrogen is so simple, it is pos- 
sible to determine the exact electronic transition rates, and there- 
fore the ratios of resulting emission lines from these transitions, as 
a function of physical c onditions in the i n terstel lar medium (see 
iMenzel & Baked dl937l) : [Baker & Menzel ( 1938 ) for the original 
theory , with updated work bv lSeatonl dl959h and lStorev & Hummed 
1 19951), and treatment s of th is t heory found in textbo o ks suc h as 
lOsterbrock & Ferlandl d2006t) or iDopita & Sutherland! d2003l) ). In 
particular, two cases exist for which the Balmer decrement has been 
determined over a range of temperatures and densities: Case A and 
Case B. Case A assumes that an ionized nebula is optically thin to 
all Lyman emission lines (i.e. lines emitted from transitions to the 
n = 1 level of hydrogen), while Case B assumes that a nebula is 
optically thick to all Lyman lines greater than Lya- (i.e. transitions 
to n = 1 from levels n = 3 and above), meaning these photons 
are absorbed and re-emitted as a combination of hya and higher 
order lines, such as the Balmer lines. These two cases will lead to 
different intrinsic ratios for the Balmer lines, with variations of the 
same order as temperature effects (for other possible "Cases" of 
emiss ion which may occur, see e.g. iFerlandll 19991 ; Luridia na et al .1 
2009). While Case B is typically assumed for determining intrinsic 
ratios, in reality the ratio in typical H n regions lies between these 
two cases, and must be determi ned using radiative transfer codes 
such as MAPPING S m (see eg lGroves et al]|2004l) or CLOUDY 
dFerlandetaill998h . 

In Table[T]we present the four strongest Balmer emission lines 
and their ratios relative to H/3 assuming Case B conditions. These 
ratios are only weakly sensitive to density, with the HayH/? ratio at 
T = 10 4 K equal to 2.86, 2.85, and 2.81 for the electron densities 
n e = 10 2 , 10 4 , and 10 6 cm 3 respectively, hence we only show the 
larger variation due t o temperature here. For a f ull ratio descrip- 
tion see Table B.7 infPopita & Sutherland! d2003l) . or Table 4.4 in 
Osterbrock & Ferland d2006l) . While these variations due to tem- 
perature and density are significant, they are still small relative to 
the effects of dust, as visible in the later sections, and hence strong 
diagnostics for the amount of reddening experienced by an emis- 
sion line galaxy. 
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2.2 Balmer absorption lines 

As discussed in the introduction, the fitting of the underlying stel- 
lar continuum is a vital step in determining line fluxes, especially 
for hydrogen (e.g. Balmer) and helium recombination lines, which 
have underlying corresponding absorption lines. The Balmer ab- 
sorption lines arise from the absorption of light by hydrogen in the 
excited n = 2 level in the photospheres of stars. The strength of the 
absorption is dependent on both the effective temperature and grav- 
ity, as they require a significant frac tion of hydrogen to be excited to 
the n = 2 level (see e.g., Table 4 in Gonzalez Delgado & Leitherer 
1999). The maximum equivalent widths occur around T e ff ~ 
9000K, corresponding to early A-type stars. 

For a simple, single-aged stellar population, the equivalent 
widths (EW) of the Balmer stellar absorption features depend 
on the age and, more weakly, on the metallicity, and vary from 
~ 2A to ~ 15 A with the maximum value occurring for stars 
aged around 500 Myr (i.e. dominated by the light from A & F 
stars). In Figure [T] we show the variation of the equivalent widths 
(EW) respectively for simple stellar populations as a function 
of age and metallicity (as labelled in the left hand of the fig- 
ure in the US diagrams). We compare the four strongest lines; 
Ho-, H/3, Hy, and US, using four models to det ermine the equiv- 
alent w idths as labelled in the upper right; the iBruzual & Charloj 
l 2003l. BC03 ) stellar population synthesis code using the MILES 
i Sanchez-Blazquez et all 120061 ; ICenarro et all 120071) and Stelib 
i Le Borgne et al ■1 120031) stellar spectral libraries, and the 2008 ver- 
sion of the Chariot & Bruzual (in prep, CB08) code with the 
MIL ES library. Also shown by the dashe d lines are the results 
from lGonzalez Delgado & Leithererld 19991) . The weak dependance 
on metallicity and the strong dependance on age, with the peak 
in EW at 0.5 - lGyr for all lines, are clearly seen for all models. 
The equivalent widths are reasonably similar between the dominant 
Bal mer lines (i.e. Ha, H/3, Hy, and Hi?) varying at most a factor of 
~ 2 dKauffmann et al.l2003l;lGonzalez Delgado & Leithererl 19991 ; 
iGonzalez Delgado et al J 19991) . 

Thus, as the equivalent widths of the stellar absorption fea- 
tures are approximately constant with wavelength while the rela- 
tive strength of the Balmer emission lines decrease rapidly with 
decreasing wavelength (i.e. for the higher order lines), stellar ab- 
sorption affects strongly the measurement of the Balmer decrement 
i.e. HafH/3, Hy/Hd This is especially so for weak emission line 
galaxies where the stellar absorption features are relatively stronger 
and only Ha is seen in emission. This relative importance of the 
effect of the stellar Balmer absorption on the emission lines is im- 
portant when considering the Balmer ratios, as discussed in later 
sections. 



3 THE SDSS SAMPLE 

Within this work we base our findings on the sp ectroscopic data 
from the seventh Data Release of the SDSS (DR7 Abazaii an et al.l 
20091). thou gh we also refer to the fourth Data Release (DR4 
Adefman-McCarfhv et alj2006 :) as well when necessary. The SDSS 
used a pair of multi-fibre spectrographs with fibres of 3" diameter. 
In most galaxies the fibres were placed as close as possible to the 
centres of the target galaxies. The flux- and wavelength-calibrated 
spectra cover the range from 3800 to 9200 A, with a resolution of 
R -1900. 

We obtain our emission line fluxes from the MPA/JHU analy- 
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Figure 1. Equivalent widths (EW) of the four dominant Balmer lines, 
Ha, H/}, Hy, and HS (as labelled), for simple stellar populations of vary- 
ing age and metallicity (increasing from top to bottom row, as marked 
on lower left in H<5 figure). Three stellar synthesis models are consid- 
ered, as indicated by th e colours in the key in the upper right; the 
IBruzual & Ch ariot (2003 , BCQ3) stellar population synthesis code using 
the MILES jSanchez-Blazquez et al.l 2006; Cena rro et aj]|2007ri and Stelib 
iLe Borgne et al J |2003|) stellar spectra libraries, and the 2008 version of 
the Chariot & Bruzual (in prep, C B08) code with the MILES library. T he 
dashed line shows the results from Gorrzalez Delgado & Leitherer J 1 999f l- 

sis of the SDSS spectroscopic samplf]. This database contains, in 
addition to the emission line fluxes, derived physical properties for 
all spectroscopically observed galaxies in the SDSS DR7. The pro- 
cedure for emission line measurement, detailed in iTremonti et al.l 
(2004), was to correct the line fluxes for stellar absorption, fitting 
a non-negative combination of stellar population synthesis mod- 
els from Chariot & Bruzual (in prep. CB08f| for the SDSS DR7 
release and IBruzual & Charlotl d2003L BC03) for the SDSS DR4 
release^]. The best-fitting stellar population model also places con- 
straints on the star formation history and metallicity of the galaxy 
(see e.g. iGallazzi et afl2 005), and has been used to estimate stellar 
masses and star-formation histories dKauffmann et all2003al) . 

The equivalent widths of the Balmer lines we use here, 
also available as part of the MPA/JHU database (listed as (line 
name)_reqw in the gal.line data file), are computed from straight 
integration over the continuum-subtracted bandpasses listed in ta- 
ble^ Note that, by definition, emission lines have negative values 
of equivalent width but for clarity in the rest of the paper we assign 
all emission lines a positive value. The continuum in this case is es- 
timated using a running median with a 200 pixel window and does 
not properly account for stellar absorption. This measurement is 
representative of the cases where a more accurate determination of 
the stellar continuum, as done with the MPA/JHU database, is not 



The data catalogues are available from 

http: //www.mpa-garching .mpg .de/SDSS/ 

2 The model spectra used were from an early version of the models and dif- 
fer from what will be eventually published. The differences from BC03 are 
primarily due to different treatment of TP-AGB stars and that the empirical 
stellar library used was the MILES library rather than STELIB. 

3 The spectra are available as part of the 
GALAXEV package, which is can be obtained from 
http: //www2 . iap. £r/users/charlot/bc2083/index.html 
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Table 2. Equivalent width bandpass 
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Line 
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Upper bound (A) 
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4101.73 


4092.0 


4111.0 


Hy 
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H/3 


4861.32 


4851.0 


4871.0 
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6562.80 


6553.0 


6573.0 



possible, such as with low resolution or low S/N data, and helps 
characterize the effects of stellar absorption on the lines. 

As we concentrate on emission line galaxies in this work, in 
particular galaxies with measurable Balmer emission lines, we have 
placed cuts on the signal-to-noise (S/N) of the Balmer emission 
lines using the uncertainties given by the MPA/JHU catalogue. As 
discussed on the website, the listed uncertainties are formal, and 
likely underestimates, thus we increase the uncertainty estimates on 
the emission lines to take into account continuum subtraction un- 
certainties by the factors listed on the web site determined by com- 
parisons of duplicate observations within the SDSS sample. Specif- 
ically for the Balmer lines, we multiply the line flux uncertainty 
estimates by a factor of 1.882. From the full SDSS emission line 
galaxy sample we define three galaxy samples depending on the 
lines and ratios being examined; SNCHasH/?), SN(Ha-,H/3,Hy), and 
SN(Hor,H/3,Hy,Ho'), where we require a S/N> 3 in two or more of 
the four strongest Balmer lines (Her, H/3, Hy, & US). The S/N cuts 
are dominated by the weakest line in each sample due to the strong 
decrease in relative flux for the higher order lines, thus the inclusion 
of each higher order line biases the samples to higher equivalent 
widths of the Ha emission line. As shown in figure [2] beginning 
from the full galaxy sample of SDSS DR7 (~ 928,000 galaxies), 
approximately half are emission line galaxies (~ 510,000, as mea- 
sured by the presence of Ho- in emission), with a broad spread of 
equivalent widths peaking at ~ 20A (as measured from the local 
continuum, not corrected for stellar absorption). As each higher 
order Balmer line is included, the sample rapidly decreases and 
is biased to higher equivalent widths, with the SN(Hci',H/3,Hy,Hd') 
sample limited to -120,000 galaxies with EW(Ha)> 10A, with 
a distribution peaking at 32A (SN(Ha-,H/3) has -392,000 galax- 
ies, while SN(Ho',Hy3,H7) has -241,000 galaxies). As the emission 
line EW(Hcf) in a galaxy spectrum can be considered a proxy for 
the specific star formation rate (the current star formation rate rel- 
ative to the total stellar mass, sSFR=SFR/M.) of a galaxy, the bias 
in EW(Ho') means a bias to more "starforming" galaxies, which 
means a bias to lower-m ass, lower-metallicity, bluer galaxies as 
show n in previous works ([Brinchmann et al. 20041; iTremonti et"al] 
|2004|) . This bias needs to be kept in mind when considering the 
diagrams and analysis in this work. 

Figure [2] also reveals the limitation of low-resolution spec- 
troscopy in finding all emission line sources, and thus the limita- 
tion of applicability of the method we explore here. Approximately 
5% of the SN(Ha) sample, and even 0.7% of the SN(Ha,H/3) sam- 
ple, actually have EW(Ho') less than zero (i.e. the emission line is 
lost in the stellar absorption feature). These sources would never be 
picked up as emission line galaxies in low-resolution spectra, and 
the use of an emission-line equivalent width Balmer decrement to 
determine the attenuation would return spurious results. 

Note that we have not included any cuts on redshift or the 
type of emission-line galaxy as in previou s works on SDSS emis- 
sion line galaxies (e.g. iKewlev et al.ll2006h . Such redshift cuts are 
necessary to make certain that aperture effects do not play a part 
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Figure 2. Distribution of the Ha emission line equivalent widths for the 
samples considered in this work. All four histograms are normalized to their 
peak value, with the total number of galaxies in each sample listed in the 
key in the upper left. The solid histogram shows all emission line galaxies 
in SDSS DR7 (defined by the presence of Ha in emission), and the three 
dashed curves show the distributions of samples defined by S/N cuts in the 
Balmer lines (as labelled in the upper left) considered in this work. 



in the derived galaxy properties (b y sampling only a s mall, biased 
part of the galaxy, as discussed in lKewlev et alj|2005l) and to pre- 
vent luminosity biases at the higher redshifts. Yet as we care only 
for the derived Balmer emission line fluxes and stellar equivalent 
widths these issues do not strongly affect our findings. 

Separating emission-line galaxies by class (i.e. star-forming or 
Active Galactic Nucleus dominated) is necessary when examining 
the structure of the forbidden emission lines (e.g. [O m]/15007A) 
which depend strong ly on the dominan t ionization mechanism 
in the gas (see e.g. Kew lev et alj l200rj) . and when comparing 
physical galaxy properti es with emission line properties (see e.g. 
Kauffman net alj l2003allbh . However, as discussed in the previ- 
ous section, the relative strength of the Balmer lines depend only 
weakly on local conditions, and will vary little in their intrin- 
sic ratios between being photoionized by AGN or by OB stars 
(i.e. Ha/Hy3 should be - 2.86 in star-forming galaxies and - 3.1 
in galaxies ionized purely by an AGN). Thus, while this difference 
is significant and will have some bearing on the work in this pa- 
per, it is secondary to the effects of dust attenuation. Only when 
collisional heating dominates the atomic gas and the excitation of 
hydrogen, such as in shocks or clouds in hot gas, can the in trinsic 
ratios be significantly different (see e.g. iFerland et al.ll2009l) . but 
these processes are not expected to dominate most galaxies within 
our sample. An additional reason to include AGN is that, in accor- 
dance with our main aim, separating out the contribution of AGN 
from lower S/N samples may prove problematic as the weak diag- 
nostic lines become lost within the noise. 

Approximately 2% of the full SDSS sample considered here 
are duplicate observations of galaxies (i.e. - 4% of the sample are 
pairs). We have not removed these from the sample so as to include 
the intrinsic scatter due to observational uncertainties in the subse- 
quent analysis. 



4 THE EQUIVALENT WIDTH BALMER DECREMENT 

The issue of simply using directly the equivalent widths of 
the Balmer lines as proxies for the line fluxes when calcu- 
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Figure 3. The variation of the equivalent width based Balmer decrement 
(log[EW(Ha)/EW(H/J)]) versus the stellar continuum-subtracted flux based 
Balmer decrement for the SDSS SN(Ha\H/?) sample. All axes are normal- 
ized by the intrinsic Balmer ratio, (Ha/HyS)i = 2.86. In the left diagram, 
the EW(H»/H/?) has not been corrected for the difference in continuum flux 
at the Ho- and H/3 wavelengths (F^fHaO/F.jfH/S)), while in the right this is 
included. The upper left hand corner shows the median uncertainties for the 
sample, and the straight line indicates a 1:1 relation. The colours and asso- 
ciated colourbars indicate the median H<5 A (, S and log[EW(Hff)] respectively 
in each pixel. The top axes give the resu lting E(B - V) from the balmer 
decrement assuming fhe lO'Donnef] jl994l) Galactic extinction curve. 



lating the Balmer decrement can be seen in figure [3] which 
shows the spread of the equivalent width based Balmer decre- 
ment (log[EW(Ha)/EW(H/3)]) against the "true" Balmer decre- 
ment determined from the stellar continuum-subtracted line fluxes 
(log(Ha/H/3)) for the SN(Ha,H/?) sample. Both the EW Balmer 
decrement and the flux Balmer decrement have been normalised 
to the intrinsic ratio of 2.86, appropriate for a low density gas of 
T = 10 4 K. 

One of the first obvious issues to be corrected for is the vari- 
ation of the underlying stellar continuum between the Ha and 
H/3 wavelengths. In the left diagram of figure [3] we show the 
distribution of log[EW(Ha)/EW(H/?)], uncorrected for the contin- 
uum flux variation, while on the right the more accurate form 
of the EW Balmer decrement is used: log[EW(Ha)/EW(HyS)] + 
log[F, t (Hci')/F / ((H/3)], where F.^Ha) is the continuum flux at Ha, 
determined from a 200 pixel median smoothing of the emission- 
line subtracted continuum. 

When uncorrected for the underlying continuum variation 
there is a clear systematic offset of the EW Balmer decrement from 
the 1 : 1 relation of ~ 0. 1 dex. Correcting for the continuum variation 
removes this offset, yet a significant spread remains. This spread is 
due to the effect of the stellar Balmer absorption features. Without 
these, EW(Ha)xF /1 (Ha) should be, by definition, the flux of the 
line. The colors in figure [3^ indicate the median H6 absorption in- 
dex (Ht^ JWorthev & Ottavianil 19971 : iKauffmann et alJ2003al) of 
the sample in each pixel. As the figure shows, while the absolute 
strength of the stellar Balmer absorption features does play a part 
in the observed offset of the SDSS galaxies' EW Balmer decre- 
ments, the dominant mechanism for the offset and spread is the 
relative strength of the stellar absorption features to the emission 
lines. This can be seen by the distribution of the equivalent width 
of the Ha emission line indicated by the colours in figure[3p, where 
there is a clear gradient of decreasing EW(Ha) with increasing off- 
set from the line. As the emission lines become weaker overall, the 
stellar absorption features, which are <10A as discussed in section 
[2] obscure a greater fraction of H/3 relative to Ha and therefore lead 
to a larger offset. 

As discussed in the introduction, the best way to compensate 
for the effect of the stellar absorption features on the emission lines 
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Figure 4. 2D histograms of the distribution of the SDSS SN(Ha,H/3) 
sample galaxies' equivalent-width based Balmer decrements 
(log[EW(Ho')/EW(H/J)]), including a constant correction for stellar 
Balmer absorption, versus the stellar continuum-subtracted flux based 
Balmer decrements. Each figure shows a different bin of Ha equivalent 
width, as labelled in lower right of each plot (in A). As in figurefj] both axes 
are normalized by the intrinsic Balmer ratio, (Ha/H/3)i = 2.86. The colours 
indicate the log of the number density in each pixel, as labelled by the 
colour bar on the right, with pixels with less than 5 galaxies excluded. The 
error bars in the upper left indicate the median uncertainty for each sample, 
with the total number of galaxies listed in the lower right of each plot. Note 
that the y-axis has a different correction for stellar Balmer absorption, R, 
for each plot as indicated in the lower left of each plot, and that galaxies 
with Hy8 still in absorption after correction (i.e. EW(Hy8) + R < 0) have been 
excluded from the sample. 



is to fit the stellar continuum as done within the MPA/JHU SDSS 
database. However, when only poor quality spectra are available 
such as for high redshift galaxies, the determination of the Balmer 
absorption features may be unreliable. One possible approach when 
faced with low resolution spectra is to assume that the absorption 
equivalent width is constant for both Balmer lines across the whole 
sample. While figure[T]clearly shows that the absorption EW is not 
the same for all the Balmer lines, it provides a first step when infor- 
mation is sparse and uncertainties large. When a constant Balmer 
absorption correction R is assumed for both EW(Ha) and EW(H/3), 
a correction factor of R = 4A is determined when the offset of the 
SDSS galaxies' EW Balmer decrements to the measured Ha/H/J 
ratios is minimized using an error-based weighting. The inclusion 
of this simple correction factor improves the situation when com- 
pared to that shown in figure [3] but with a still significant scatter 
of a ~ 0.1 dex around the expected value and an extended tail 
of objects towards lower values. Both the scatter and the tail arise 
due to the assumption of a constant offset (i.e. Balmer absorption) 
for the whole sample. The value determined is biased towards high 
EW(Ha) galaxies, as these galaxies both dominate the sample and 
have lower uncertainties, as discussed in section[3] 

When split into bins of different EW(Ha), more accurate fits 
with differing correction factors are obtained. In figure|4] we show 
the fits for the galaxies split into four bins; 0.5 < log(EW(Ha)) < 
1.0, 1.0 < log(EW(Ha)) < 1.5, 1.5 < log(EW(Ha)) < 2.0, and 
2.0 < log(EW(Ha)) < 2.5. The number of galaxies in each bin is 
listed in the lower right of each figure, and the median uncertainties 
are shown by the error bars in the upper left. Note that galaxies with 
H/J still in absorption after the correction factor R is applied have 
been excluded, but that this is less than 0.2% of the sample in each 
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bin. As can be seen in the figure, the median uncertainties increase 
quickly with decreasing emission line equivalent width. It is for 
this reason that galaxies with log(EW(Ha)) < 0.5 have not been 
included here. 

The y-axis for all four plots includes a correction for stellar ab- 
sorption; log([EW(Hor) +i?]/[EW(H/3) +R]) +log[F A (Ha)/F, ,(H/3)]. 
As for the full sample, we determine the correction factor, R, for 
each binned sample by finding the value that leads to the minimum 
offset from the 1:1 relation with the uncertainties giving lcr off- 
sets from this relation. The values determined are R = 3.1, 3.5, 
4.1, and 4.1A for each increasing bin of EW respectively, as indi- 
cated in the lower-left of each plot in figure[4] The lcr uncertainty 
around R for each bin is approximately 1.0 (slightly less for the 
EW(Hci-)> 100A bin). The lowest bin has an uncertainty of 1.5, 
though the probability distribution for R is slightly skewed to higher 
values, arising from the offset visible in|4^, discussed below. The 
Balmer decrement determined from the EWs is significantly better 
than for the full sample, and especially so when no correction is 
included (figure[3]>, with dispersions of lcr ~0. 11, 0.06, 0.05, and 
0.04 dex around the 1 : 1 line respectively for the 4 binned samples. 

The lowest EW(Hcr) (figure [4^) sample appears by eye to be 
slightly offset form the line. A reasonable hypothesis is that this 
offset arises due to our simple assumption of a constant correction 
R to both EW(Hcr) and EW(H/3), whereas it is clear from[T]that the 
absorption EW(Ha) is typically less than the absorption EW of H/? 
by approximately a factor of 0.6 on average. However, changing 
the correction factor of EW(Ha) to 0.6R and redoing the fit does 
not remove this offset. On closer examination it is clear that this 
offset is due to a biasing of the fit to the high signal-to-noise data, 
which predominantly occur at high values of the Ha/H/3 ratio due to 
measurement biases (i.e. there is a clear gradient in EW(Hci-) from 
top to bottom in figure [4^). Assuming uniform weighting for the 
fit (i.e. ignoring the errors in EWs) gives a value of R = 3.0, well 
within the large uncertainties for R. For the other figures, assuming 
an offset correction factor for EW(Ho') of 0.6R, results in R = 3.5, 
3.7, and 3.7A in terms of increasing EWs, with similar dispersion 
around the relations. The results are within the uncertainties for R 
when assuming a constant correction, but in all cases indicate the 
necessity for the correction of stellar absorption to the emission 
lines of a factor of ~ 4A. 

For figures gja, b, and d (0.5 < log(EW(Ha)) < 1.0, 1.0 < 
log(EW(Ha)) < 1.5, and 2.0 < log(EW(Hff)) < 2.5 respectively) 
the observed scatter is less than the median uncertainty. Only for 
figure [4J; does there appear to be a significant, low number scat- 
ter above the line (note the log scale density in figure|4](. However, 
when examined, the scatter in this diagram, and also in figures [4^, 
b, and d, is correlated with the uncertainty in the determined Her 
and H/? lines, with the median scatter in each pixel increasing sig- 
nificantly the further from the line. Thus the median uncertainty for 
the outliers is greater than the median uncertainty of the sample as 
a whole in figured- 
While figure[4]demonstrates that it is possible to determine the 
Balmer decrement to some accuracy from emission line equivalent 
widths, the determination of the correction is dependent upon the 
measurement of three quantities; EW(Ha), EW(H/3), and the flux 
ratio, F / i(Ho')/F, 1 (H/3). While the former two will be observable in 
strong emission line galaxies at high redshift, the flux ratio may 
prove problematic to measure from spectra. However this ratio is 
closely tied with the observed optical colours of the galaxy. In fig- 
ure|5]we show the distribution of the continuum fluxes measured at 
the Her and H/? wavelengths (F, 1 (Ho')/F^(H/3)) against the restframe 
g - ;• colour as measured from the SDSS fibre spectrum within the 
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Figure 5. Distribution of the full SN(Ho-,H/3) SDSS sample ratio of con- 
tinuum fluxes measured at the Ho- and H/S wavelengths (F ,i(HoO /F ,i(H/3)) 
against the restframe g—r colour as measured from the SDSS fibre spectrum. 
The solid line with diamonds shows the median F,((Ho')/F^(H J 6), while the 
dashed lines show the lcr dispersion. Median uncertainties are indicated by 
the eiTor bars in the upper left. The thick dot-dot-dashed line indicates the 
best fit linear relation, with y = -0.26 + 0.39(g - r) . 

full SN(Ha,Hy8) sample. A linear fit to the correlation in this figure 
returns 

log (F,,(HqO/F,(H/?)) = -0.26 + 0.39(g - r) 6pec , (1) 

with a standard deviation of tr ~ 0.015 dex around this relation. We 
use the rest-frame g — r colour as a proxy for stellar continuum, but 
other colours such as r - i provide similar constraints. In the case of 
low-S/N spectra, the rest-frame colour could be from SED fitting 
to broad-band magnitudes. 

Thus, in summary, the Balmer decrement for a low-resolution, 
strong-emission-line spectrum can be measured from the emission 
line equivalent widths and colours alone with; 

/ EW(Hcr) + 4 1 \ 

(2) 

with a scatter around this of cr ~ 0.05 dex, or ~ 0.3 m ag in 
A v , assuming a Galactic extinction law (e.g. lO'Donnellll 19941) . For 
weaker emission line galaxies ((i.e. EW(Hcr)< 30A), a lower offset 
(R ~ 3.5) should be used, as shown in figure|4] However, given the 
larger uncertainties and greater dispersion seen at lower EW(Hor), 
a correction factor of R ~ 4 can be used for the full sample with a 
0.1 dex scatter (~ 0.7 mag in A v ) and an extension to lower values 
(i.e EW Balmer decrement underestimate) due to low EW systems. 

One final note on this relation: as seen in figure [3] there is a 
strong bias in the sample of Balmer decrement with other galaxy 
pro perties, as discussed in detail in several othe r SDSS papers (see 
e.g. lKauffmann et al.l2003al : lGarn & Bes^l2010h . Thus, the relation 
shown above includes a combination of both galaxy type as well 
as variation in extinction. The only way to remove fully this effect 
is to match pairs of galax ies in as many pro perty types excluding 
extinction, such as done in lWild et al. I J201 lbl) . Unfortunately when 
applied to the sample here, it was found that the range in extinction 
was not large enough to properly determine the relation. However, 
even given these uncertainties, this relation should still hold at sev- 
eral redshifts as high attenuations are on average associated with 
high gas masses, and thus high star formation rates and similar un- 
derlying continua at all redshifts. 
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Figure 6. The distribution of Hy/H/3 ratios versus HafHfi ratios for the 
SDSS SN(Ha,H/?,Hy) sample, with the number density of each pixel in- 
dicated by the bar on the right and pixels with less than 10 galaxies not 
shown. Both axes are normalised to their respective Case B ratios with 
(Ha/H8)i = 2.86 and (Hy/H/3)i = 0.468. The median uncertainty of the 
sample is indicated by the error bars on the right. Overplotte d are the vari- 
ation expected from three different attenuation laws; Chariot & Fall (2000, 
CF00, diamonds), lO'Donnell jl994 0'D94, triangles), an JCalzetti et alj 
1 2000, CalOO, plus signs) (see text). Each symbol represents a step of 0.5 in 
Ay, up to A y = 3. Inset: A zoomed in version of the figure, showing the po- 
sition of the zero point when the intrinsic ratios ((Ha/H/3)i and (Hy/H/3)i) 
are taken to be at T=5,000K (diamond: 3.04, 0.458), at 10,000K (circle: 
2.86, 0.468), and at 20,000K (tri angle: 2.75, 0.475). T he solid line shows 
the variation of the ratios with the Calzetti et al. 1 2000) attenuation law. 



5 STELLAR ABSORPTION EFFECTS ON THE 
EMISSION LINES 

One issue with the previous section is that we assume throughout 
that the stellar-continuum subtracted emission line fluxes within the 
MPA/JHU database are correct. While the overall fits to the stellar 
continuum are impressively good with a median^ 2 = l.Olperpixel 
across the sample, there are appear to be remaining issues around 
the Balmer lines. In the following we concentrate on the SDSS DR7 
Balmer emission-line fluxes corrected for the underlying stellar ab- 
sorption features from the MPA/JHU database, and explore their 
uncertainties using the known intrinsic values and commonly used 
attenuation and extinction laws. 



5.1 The issue with H/3 

When considered alone, the ratio of Hcr/H/3 cannot indicate prob- 
lems with the measurement of the lines involved unless it is signif- 
icantly below the expected value of the unattenuated ratio. This is 
because the larger values of the emission line ratio can be caused by 
attenuation by intervening dust, with the intrinsic ratio dependent 
the emitting gas density and temperature (as discussed in section 
O. However by examining several of the Balmer lines at once these 
dependencies can be accounted for. 

Figure [6] shows the variation of the Hy/H/? ratio against the 
HafHfl ratio for the SN(Ha,Hj8,Hy) SDSS sample. Both ratios have 
been normalized to their Case B, T=10,000K, n e = 100 cm~ 3 in- 
trinsic ratios ((Ha/H/% = 2.86, (ny/H/3\ = 0.468). Immediately 
obvious in this figure is the offset of the sample from the zero point, 
indicating most SDSS galaxies undergo some attenuation (as seen 



in the previous plots), with the correlation between the two ratios 
as expected from the reddening laws applied to the intrinsic ratio. 

The three different lines overplotted show the effect of three 
different attenuation laws commonly assumed in the analysis of 
galaxies. For all three lines, the symbols indicates steps of 0.5 in 
Ay, up to A v = 3. 

ThelO'Donnelll dl994) law (OD94) is an updated version of 
the ICardelli et alj dl989h fit to the average extinction law in the 
Galaxy, thus least representative of the integrated emission from 
the SDSS galaxies, which will suffer attenuation due to the mixture 
of emitting source s and absorbing medium. However, as discussed 
in iKennicutt et alj d2009h and can be seen in figure [6] the use of 
a foreground dust screen with galactic extinction is indistinguish- 
able from the other laws, especially given the uncertainty within 
the SDSS sample. We assume a total to selective V-band extinction 
of R v = 3 .1, the average value in our galaxy. 

The ICalzetti et al] fcOOd) attenuation law (CalOO) was ob- 
tained from the continuum and Balmer decrement of local actively 
star-forming galaxies, thus matching the high EW(Ha) galaxies in 
the sample. Note that as only ratios are analysed here, the differ- 
ence between the colour excess (E(B - V)) of the stellar continuum 
and nebular lines noted by Calzetti et al. is effectively sc aled out. 
The R v used here is 4.05, as given bv lCalzetti etal] ( 2000) from the 
comparison of the observed infrared flux to that predicted from the 
obscurati on of the optica l -ultra violet light. 

The ICharlot & Falll J200ch attenuation law (CFOO) is a more 
simple, empirical law put forward to all ow for the different colour 
excesses and attenuations observed by ICalzetti et al. | j2000h be- 
tween the nebular emission lines and stellar continuum. It breaks 
the attenuation into two components; the 'diffuse ISM' component 
that describes the effective obscuration of all stars in a galaxy by 
the diffuse dust, and the 'birth cloud' component that describes the 
additional extinction suffered by the Hn regions from which the 
nebular emission lines arise, giving 



Ay 



(3) 



where A v = 5500A. Th e exponent of -0.7 fo r the diffuse ISM was 
empirically derived by ICharlot & Fall j2000) with a comparison 
of nearby galaxies, while the -1.3 exponent for the birth clouds 
was chosen to match the extinction within our own Galaxy. The 
parameter /j indicates the fraction of the attenuation suffered by the 
ne bular lines by each com ponent. We assume n = 0.3 here, as used 
bv lWildetaU d2007l) and IWild et at] feOllj) in their analyses of 
SDSS galaxies. 

There are three things of note to take from figure[6] One, given 
enough precision in the data, the SDSS galaxies should be able to 
distinguish between these three different attenuation laws, and pro- 
vide an answer on which law is best (or least bad) to apply to an 
ensemble of galaxies. Such work has been done before f or com- 
parin g galactic extinction laws using planetary nebulae (Phillips 
120071) . Note that the different R v between the CalOO and OD94 
laws is what causes the difference in expected Ha/HfJ for the same 
A v , while the CFOO and CalOO have similar A v as determined from 
Hff/H/3, but not from Hy/H/3. Two, that the scatter of the SDSS 
galaxies is large around the three laws, preventing this possibility, 
though this scatter is not significant when compared to the median 
uncertainty as shown by the error bars. Third, and most importantly, 
while the slope of SDSS galaxies matches that given by the atten- 
uation laws, there is a systematic offset of the SDSS galaxies when 
compared to all three attenuation laws. This offset is significant, 
and cannot be explained by assuming more extreme (and therefore 
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Figure 7. The same plot as figure[6] except using the SDSS DR4 MPA/JHU 
catalogue. The same SN cuts on Ha, H/3, and Hy are applied, with ~ 160,000 
galaxies in the sample. 



less likely) attenuation laws or by assuming large values of R v , as 
the zero point itself appears to be offset. 

Neither can this situation be remedied by assuming different 
values for the unattenuated Balmer ratios. As discussed in section[2] 
the intrinsic Balmer emission-line ratios are sensitive to the temper- 
ature of the ionized gas from which they arise, and, more weakly, 
to the density of the gas as well The inset in figure [6] shows a close 
up of the zero point of figure [6] i.e. galaxies with little or no at- 
tenuation. Over this are plotted 3 symbols indicating the position 
where the zero point would be for 3 different average Hn region 
temperatures; 5,000K, 10,000K (assumed within the figure), and 
20,000K. All assume Case B ratios, and a typical H n gas density 
of n c = 10 2 cm~ 3 . These 3 temperatures encompass the range of 
temperatures expected, with typical solar metallicity Hn regions 
having T c ~8,000K. What this figure demonstrates is that the vari- 
ation in intrinsic Balmer ratios is small relative to both the effects of 
attenuation and the observed offset, and that the variation in intrin- 
sic ratios is in the same sense as that du e to the effects of attenuation 
(as shown by the ICalzetti et al. (2000) law), thus the intrinsic ratio 
cannot be the cause of the offset. 

A possible cause for this offset can be found when the 
MPA/JHU DR4 catalogue is examined instead. Using the same SN 
cuts on He, HyS, and Hy, figure [7] shows the same plot as figure [6] 
for the DR4 sample. In most ways this figure is exactly the same 
(as it should be), except in two respects: the DR4 SN(Ha,Hy8,Hy) 
sample only has ~ 160,000 galaxies in the sample compared to the 
~ 240,000 galaxies in the DR7 sample, visible in both the number 
density (colours) and in the scatter, and the DR4 sample does not 
have the offset with respect to the attenuation laws seen in figure[6] 

As mentioned in section [3] the major difference between the 
DR4 and DR7 line fluxes from the MPA/JHU catalogues is the ver- 
sion of the GALAXEV models used for the continuum fits. In some 
respects the difference seen in the figures is surprising, as the me- 
dian^ 2 of the fits to the continuum in DR7 is reduced compared to 
the DR4 fits, from^ 2 = 1.5 for DR4 to^ 2 = 1 for DR7, suggesting 
a much better fit to the spectra. 

The issue most likely lies within the H/? line region, as an 
offset in the H/3 line flux would also explain the increasing off- 
set at higher (Ha/H/J) between the SDSS galaxies and the atten- 
uation laws, and this region has been observed to be mismatched 
between models and the spectra of some globular clusters (see e.g. 
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Figure 8. a) The s ame diagram as Fig ure|6] but rotated to A y axes using the 
attenuation law of Calzetti et al. ( 200y). The x-axis is the Ay as measured 
from the (Ha/H/J) ratio, while the y-axis is the offset from the expected 
(Hy/H/?) ratio based on this Ay. b)The same as a), but the H/J line flux has 
been corrected for under subtraction of the stellar continuum by 0.35A. 



IWalcher et aI1l2009l ; IPooie et alJl2010T) . This difference in slope is 
more clearly seen in figure [8^ where we have rotate d and scaled 
figure [6] to the Ay plane using the attenuation law of ICalzetti et al] 
(2000), where the offset from the x-axis is more clearly seen. Simi- 
lar results are seen if another attenuation law is used. In this frame, 
the y-axis is the offset from the expected (Hy/H/?) ratio based on 
the Ay as determined by the (Ha/H/?) ratio. 

By "correcting" for the incorrectly subtracted stellar H/J equiv- 
alent width we can fix figure [8^. This is what we have done in fig- 
ure[8]j, where we have added 0.35A to the emission line equivalent 
width (i.e. "True" 1^3=^8+0.35?^^)). The correction of 0.35 A 
was determined by minimizing the offset from the Calzetti law (this 
value depends only weakly on the choice of attenuation law). This 
means that the stellar absorption equivalent width of H/? is system- 
atically underestimated by 0.35A in the CB08 continuum fits to 
DR7 SDSS spectra. 

The fact that this 0.35A underestimation is systematic is in- 
teresting, as we would expect that any error would depend on the 
strength of the Balmer features, either measured through the H5 Abs 
or the emission line equivalent widths, both of which do correlate 
with Ha/H/? as seen in figure [3] , yet no correlation is observed. 
What the underlying cause of this systematic underestimation is 
unknown, yet it must be taken into consideration when determin- 



ing the A v from the Balmer decrement in the SDSS DR7 MPA/JHU 
catalogue. It is possible that this issue arises due to the misclas- 
sification of the spectral resolution o f the MILES library (as dis- 
cussed in lFalcon-Barroso et alJuOllh in the implementation in the 
CB08 code, but this issue is now known and currently under in- 
vestigation. This investigation goes beyond the scope of the work 
presented here but we note that the models used here are early 
versions of the CB08 library and these issues are expected to be 
solved within the to-be-published models. When the new Ay is cal- 
culated from the Htt/ H/3 ratio in cluding the systematic 0.35A offset 
(using e.g. Calzetti e t alJl200oL law), a mean difference of -0.07 
magnitudes is found with the uncorrected A v estimates, increas- 
ing slightly at high er Ay. This s uggests that previous DR7 Ay es- 
timates, such as in lGarn & Best! l l2010l) . are overestimated by this 
value. 

Similarly this 0.35A correction must be included in our EW 
Balmer decrement (equation[2]l leading to a new equation, 

logdfa/H/J) = log( EW ( (H/? ; + 44 ) + (-0.26 + 0.39(g - r )rest ) , 

(4) 

which more closely matches our expectations of different stellar ab- 
sorption equivalent widths between Ho- and H/3. For weaker emis- 
sion line galaxies, the offset would be smaller than 4. 1 , as shown in 
figure|4]and discussed at the end of section|4] 



The Balmer decrement of SDSS galaxies 9 




0.0 0.5 1.0 1.5 2.0 

A, (Hy/Ha) (CalOO) 



Figure 9. The distribution of H<S/Ha versus Hy/Ha for the 
SN(Hff,H/?,Hy,H,S) sample of SDSS DR7 galaxies (-121,00 objects). 
As in figure[8] the fi gure i s rotat ed and scaled to the Ay plane as determined 
from the Calzetti et al. 12000) law, indicated by the straight line. The 
number of galaxies in each pixel is indicated by colours (as labelled), with 
pixels with less than 5 galaxies not shown. 



5.2 The issue with US 

While examining the issue in H/3, a similar issue was found for HS 
that we present here as a curiosity. When HS/Ha versus Hy/Hais 
plotted, using the SN(Ha,Hfi,Hy,HS) sample of SDSS galaxies and 
avoiding the problematic H/3 line, the tight correlation between 
these two ratios, matching closely that expected from the attenu- 
ation by dust. However, upon closer examination, a systematic off- 
set is observed between the galaxies and attenuation laws. As in 
the previous subsection on H/3, the offset is clearer when rotated to 
the A v plane, which is what we show in figure [9] where the x-axis 
is the A v determined from Hy/Hor using the ICalzetti et al.l bOOd) 
law, with the v-axis the offset from this A v when determined from 

the HS/Ha ratio. The median offset is 0.05 for most of the A v 

range shown here, meaning that the Ay determined from HS/Ha is 
lower than that determined from Hy/Hor. The offset is seen for all 
attenuation laws considered here. 

More importantly, the offset is also observed in the DR4 sam- 
ple, though with greater uncertainty due to low number statistics. 
As with the Hf} line there appears to be no correlation of the offset 
with emission line equivalent widths, or stellar age determin ants 
like H<5 Abs or the D„4000 index (see iKauffmann et al j|2003al , for 
definitions of these indices). Neither does it appear to be correlated 
with Ha/Hp or other emission line or attenuation tracers. Thus, due 
to the lack of difference between DR4 and DR7, the strong EW 
bias of the sample (as shown in figure O, and the fact that the un- 
certainty is dominated by the H<5 line, it is still not known what 
exactly causes this offset. It is most likely an issue due to the under- 
lying continuum, but an investigation into the stellar models goes 
beyond the scope of this work. Thus we present this issue for now 
as a curiosity and a cautionary note of the level of systematic un- 
certainties in determining weak line fluxes from the SDSS sample. 



6 CONCLUSION 

We have examined the possibility of using equivalent widths of 
the Balmer emission lines to determine the Balmer decrement, and 
hence attenuation, of a galaxy. Using the Sloan Digital Sky Sur- 
vey we were able to determine a statistically representative relation 
between the continuum-subtracted Balmer emission line flux ra- 
tio and the equivalent widths (EW) of the Balmer emission-lines 
combined with a rest-frame colour, correcting for the effects of the 
stellar absorption features: 

/EW(Ha) + 4.1 \ 
l°g(H«/H0) = log ( EW \ H J ) + 4A j + (-0.26 + 0.39(g - r)) , (5) 

for galaxies with EW(Ha)> 30A, with a scatter of lcr ~ 0.06 dex, 
or 0.4 mag in Ay, indicating the possible variation for individual 
objects. For galaxies with EW(Ho-)< 30A smaller correction fac- 
tors (3.5 for EW(Hor), 3.8 for EW(H/3)) should be used. However, 
given the scatter at low EW values, the equation above can be used 
above for all galaxies allowing for a much greater uncertainty in 
the final Balmer ratio or Ay determined. 

In addition, by comparing the Balmer decrement (Ha/H/3 
versus Hy/H/3) we discovered that the H/3 emission line equiva- 
lent width (and hence flux) is underestimated by 0.35A in the 
JHU/MPA DR7 SDSS database, due to an issue in the H/3 region 
of the 2008 version of the Chariot & Bruzual stellar population 
synthesis code GALEXEV. This leads to an overestimation of the 
atten uation of the SDSS g alaxies of 0.07 magnitudes in Ay assum- 
ing a lCalzetti et al.l j200Ct) attenuation law. 

Finally, we also discovered a strange offset in the HS emis- 
sion line fluxes observable in both the DR4 and DR7 releases of 
the MPA/JHU database which we present both as a curiosity and as 
a warning on the underlying issues in interpreting weak-line emis- 
sion lines in a statistical sample 
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